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In  the  Uniform  Segment  Method  the  spatial  domain  is  divided  into  ajfinite  number 
of  prismatic  sections  within  which  the  Euler  equation  is  appl ied.-^ Solutions  to 
this  equation  yield  functional  relationships  for  mode  shapes  in  the  form  of 
trigonometric  and  hyperbolic  functions.  The  overall  structure  is  modelled  as  a 
free-free  beam  with  intersegment  continuity  assured  by  the  matching  of  displace¬ 
ment,  slope,  moment,  and  shear  transfer  at  segment  boundaries.  Applying  the 
free-free  boundary  and  continuity  conditions  results  in  a  system  of  equations  in 
the  unknown  frequencies  and  coefficients.  The  solution  of  this  system  yields  thej 
natural  frequencies  and  mode  shape  coefficients  within  an  arbitrary  constant. 

The  main  advantage  of  using  this  method  over  finite  elements  lies  in  the  reduced 
number  of  degrees  of  freedom  needed  to  model  the  structure.  Real  prismatic 
sections  model  one  to  one.  This  is  not  the  case  in  finite  element  methods. 
Additionally,  a  number  of  the  transient  loads  known  to  drive  gun  vibrations  are 
functions  of  the  mode  shape  derivatives.  Since  these  functions  are  represented 
analytically,  which  is  a  characteristic  of  the  Uniform  Segment  Method,  the  need 
to  calculate  derivatives  by  numerical  algorithms  is  eliminated.  Thus,  exactness 
in  derivative  calculations  is  assured. 
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BACKGROUND 


At  the  Fourth  U.S.  Army  Symposium  on  Gun  Dynamics,  which  convened  in  May 
1985,  the  author,  in  collaboration  with  H.  J.  Sneck,  submitted  a  paper  (ref  1) 
in  which  a  method  for  modelling  the  flexural  vibrations  of  a  tank  gun  using  the 
normal  modes  technique  was  proposed  and  developed.  In  this  report,  a  number  of 
shortcomings  in  both  the  model  and  analysis  techniques  are  cited  and  will  be 
reviewed  briefly. 

The  first  and  most  striking  model  shortcoming  lies  in  the  use  of  an  axially 
prismatic  beam  as  representative  of  the  gun  tube  knowing  full  well  that  a  gun 
possessing  a  prismatic  barrel  does  not  exist.  The  reason  for  the  choice, 
however,  lies  in  the  fact  that  although  guns  are  axially  non-uniform,  their 
deviation  from  the  prismatic  condition  is  small  when  compared  to  their  overall 
length.  For  example,  when  viewed  from  a  distance,  any  gun  appears  to  be  of  uni¬ 
form  cross  section.  A  less  important  reason  for  this  choice  was  the  ease  in 
which  this  type  of  model  could  be  developed.  For  the  prismatic  case,  a  single 
four-term  mode  shape  equation  results  for  each  natural  frequency  considered. 

The  second  involves  the  choice  of  support  conditions.  In  this  initial 
effort,  the  supports  were  modelled  as  linear,  bi-directional  springs  rigidly 
attached  to  the  breech  end  (origin  of  axial  coordinate)  of  the  beam.  One  spring 
resisted  the  transverse  displacement,  while  the  second  restrained  the  rotational 
motion  of  the  breech  end.  This  type  of  support  model  is  not  exact,  but  rather 
an  approximation  of  what  actually  exists  on  fielded  weapons.  In  addition,  sup¬ 
ports  possessing  no  clearance  are  not  feasible  from  an  assembly  or  operational 
standpoint  and  those  possessing  linear  response  properties  regarding  force/ 
deflection  are  not  likely. 


References  are  listed  at  the  end  of  this  report. 


The  third  model  shortcoming  involves  the  use  of  the  Eu^er-Berncul  1  ->  beam 
equation  without  a  strong  claim  regarding  its  applicability  for  gun/beam  type 


structures.  It  was  chosen  based  upon  the  results  of  an  order  of  magnitude  study 
on  the  equation's  coefficients  for  the  model  of  a  gun  system  similar  to  the  M60. 
A  more  viable  approach  comparing  the  results  of  a  modelling  effort  using  alter¬ 
native  equations  would  be  a  better  means  of  justification. 

A  fourth  condition  involves  the  identified  excitation  sources  (driving 
loads)  developed  into  the  model.  There  were  five  sources  which  included  recoil 
inertia,  pressure  curvature,  projectile  trajectory,  projectile  rotational  imbal¬ 
ance,  and  eccentricity.  These,  by  far,  are  the  most  important  loading  sources, 
however,  a  more  versatile  model  employing  additional  sources  would  be  more 
realistic. 

All  of  the  above  areas  are  addressed  in  the  current  analysis.  This  report 
highlights  the  detailed  improvements  which  have  been  developed  into  the  current 
computer  model  dedicated  to  gun  vibration  analysis.  The  routines  are  available 
and  running  on  Benet  Laboratories  VM/SP  Operating  System. 

FORMULATION  OF  THE  CANNON  BEAM  DYNAMICS  MODEL 
General  Modelling  Equation  and  Solution  Proposal 

Beam  dynamics  as  applied  to  gun  tubes  is  represented  in  Figure  1  where  the 
beam  structure  shown  possesses  axially  varying  geometric  properties  inferring 
that  both  bending  and  inertial  resistance  are  functions  of  position.  The  trans¬ 
verse  cross  section  is  cylindrical  and  axisymmetric,  and  in  its  static  state, 
the  axis  assumes  in-plane  curvature  due  to  the  beam's  weight,  environmental  con¬ 
dition,  and  the  result  of  manufacturing  processes.  The  features  drawn  in 


phantom  represent  the  additional  mass  of  the  breech,  bore  evacuator,  and  muzzle 
brake.  Four  types  of  forces  and/or  moments  are  represented  and  appropriately 
placed  upon  the  structure. 


CANNON/BEAM  LOADING  SCHEMATIC 

TYPE  3  TYPE  2  TYPE  4 


Figure  1.  Cannon/beam  loading  schematic. 

The  equation  which  completely  describes  the  transverse  motion  of  this  beam 
system  according  to  Thomson  (ref  2)  and  others  (refs  3-6),  is: 


N 
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where 


E  *  Young's  Modulus  of  material 


I  *  transverse  moment  of  inertia 


J  =  pitch  moment  of  inertia 


G  *  modulus  of  elasticity  in  shear 
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A  =  area  of  beam's  cross  section 
k  =  shape  factor  of  cross  section 
w  =  weight  per  unit  length  of  beam 
p-j  *  i-th  forcing  function  (total  of  N) 
g  =  gravity  acceleration 

x  =  axial  coordinate  (independent  space  variable) 
t  »  time  (independent  time  variable) 
y  =  transverse  displacement  (dependent  space  variable) 

•  =  time  derivative 

'  =  space  derivative 

This  equation  is  commonly  referred  to  as  the  Timoshenko  equation.  For  a  pris¬ 
matic  beam,  the  second  differential  operation  on  the  first  term  could  be  placed 
within  the  brackets  since  the  bending  resistance  becomes  constant.  The  second 
term  deals  with  the  translational  inertia  of  a  beam  segment,  while  the  remaining 
two  terms  model  rotatory  inertia  and  shear  deformation  effects.  All  terms  are 
on  a  per  length  basis.  The  right  side  of  the  equation  contains  the  represent¬ 
ative  driving  and  reaction  loads. 

Reduced  forms  of  this  equation  are  used  when  the  geometry  causes  some  of 
the  coefficients  to  be  relatively  insignificant.  Other  special  cases  arise  for 
the  analysis  of  prismatic  beams  for  which  all  geometric  properties  are  constant 
in  space  rendering  all  dependent  variable  coefficients  constant.  This  allows 
for  a  closed-form  solution  to  the  homogeneous  equation  as  will  be  addressed 
later.  For  'thin'  structures,  terms  involving  rotatory  inertia  become  small, 
thus  modifying  the  coefficient  of  the  third  term  and  completely  eliminating  the 
fourth  term.  For  'long'  structures,  bending  dominates  shear,  thus  both  terms 
three  and  four  may  be  eliminated  without  loss  of  modelling  accuracy.  In  this 
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simpler  form,  the  Timoshenko  equation  reduces  to  the  Euler-Bernoull i  equation, 
which  upon  the  stipulation  of  constant  cross-sectional  properties,  becomes 


N 

Ely""  +  -  y  =  5  [Pi (x,t,y,y' ,y",y' )]  -  w  (2) 

9  . 

i=l 

The  solution  technique  for  the  proposed  model  in  its  full  or  any  of  the 
truncated  forms  is  the  same.  The  homogeneous  equation  will  be  solved  by  the 
separation  of  variables  technique  yielding  the  normal  modes  of  vibration  of  a 
free-free  beam  structure.  The  model  equation  will  be  reassembled  using  the  mode 
shape  vectors,  and  upon  invocation  of  the  appropriate  orthogonality  principle,  a 
system  of  ordinary  differential  equations  (O.D.E.'s)  containing  time  varying 
amplitude  coefficients  will  result.  There  will  be  one  equation  for  each  mode 
shape  considered.  Due  to  the  nature  of  the  loadings,  this  system  of  O.D.E.'s  is 
coupled;  therefore,  numerical  matrix  procedures  will  be  needed  to  arrive  at  the 
total  solution. 

Assessing  the  Significance  of  Shear  and  Rotatory  Inertia 

The  Timoshenko  equation  represents  the  most  complete  form  of  beam  analysis 
available.  For  certain  types  of  beam  structures,  the  combined  effect  of  shear 
deflection  and  rotatory  inertia  contributes  little  to  the  results  while  their 
inclusion  is  costly  in  terms  of  solution  complexity.  Since  this  modelling 
effort  is  to  be  applied  to  a  certain  class  of  structures,  namely  large  caliber 
cannon  tubes,  these  effects  may  be  judiciously  neglected  in  favor  of  the  less 
complex  Euler-Bernoul 1 i  model.  The  reason  for  this  will  be  developed  in  the 
following  comparative  analysis  of  a  uniform  'thin'  beam. 

Consider  the  free  vibrations  of  an  unrestrained  prismatic  beam.  The 
Timoshenko  equation  (1)  becomes 


(3) 
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Ely""  +  -  y  -  (J  + - )y"  +  ( - )y  =  0 
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Following  the  method  used  by  Bozich  (ref  4)  and  identifying  the  parameters 


suggested  by  Kruzlewski  (ref  7),  this  equation  may  be  reformulated  in  a  dimen¬ 


sionless  form. 


1  d4y  d2y 

...  (...)  +  (KR2+KS2) - (1  -  KB2KR2Ks2)y  =  0 

kB*  dz4  dz2 


where 


Kg2  =  u2L4  --- 

9ei 


1  EI 


Ks2  = - 

L2  kAG 


1  Jg 


Kr2  =  --  - 


with  u  being  a  vibration  frequency.  These  three  coefficients  represent  the 


contributions  of  bending,  shear,  and  rotatory  inertia,  respectively.  Equation 


(4)  represents  the  non-dimensional  form  of  the  Timoshenko  model  which  reduces  to 


the  Euler-Bernoul 1 i  version  by  simply  setting  Ks  and  KR  to  zero. 


In  terms  of  the  non-dimensional  coordinate  z,  the  analytical  solution  to 


the  equation  is 


y(x,t)  =  (Cj  cos  az  +  C2  sin  oz  ♦  C3  cosh  0z  +  C4  sinh  /Sz) cos  wt  (6) 


where 


4  (KR2+KS2)  +  /(KR2-KS2)  ♦  4 /Kg 2 


/-  (Kr8+KS2)  ♦  /(KR2-KS2)  +  4/Kg* 


These  equations  are  valid  when  the  terms  within  the  outer  radical  remain  posi¬ 
tive.  By  applying  free-free  boundary  conditions  and  normalizing  the  mode  shapes 
with  respect  to  Cj,  the  following  characteristic  equation  for  the  natural 
frequencies  may  be  derived: 

/3 2 -a 2 

2(l-cos  a  cosh  /S)  +  ( - )  sin  a  sinh  fi  =  0  (8) 

a/3 

This  equation,  which  is  essentially  a  function  of  u>,  has  an  infinite  number  of 
roots.  It  may  be  solved  to  any  desired  accuracy  for  a  finite  number  of  them  by 
standard  root  finding  techniques.  A  set  of  computer  routines  for  solving  this 
uniform  beam  frequency  equation  was  written  for  the  assessment  analysis.  The 
effect  of  shear  deflection  and  rotatory  inertia  may  be  included  or  neglected  in 
order  that  an  evaluation  of  their  effects  may  be  made.  The  equation  coef¬ 
ficients  evaluated  for  each  frequency  are 


Cl  -  1 

/S  cos  a  -  cosh  0 

(9a) 

a  sinh  p  -  (/3/a)sin  a 

a 

(9b) 

c3  -  <-P 

a 

(9c) 

c4  =  m»c2 

p 

(9d) 

These  coefficients  and  their  respective  natural  frequencies  are  used  to  evaluate 
and  plot  mode  shapes  as  a  function  of  the  dimensionless  coordinate  z.  The 
results,  as  calculated  by  the  Timoshenko  and  Euler-Bernoull i  models  for  a  given 
structure,  may  then  be  compared  to  determine  whether  the  shear  deflection  and 
rotatory  inertia  effects  are  significant  for  the  class  of  beams  to  be  analyzed. 


In  order  to  assess  the  effect  of  shear  and  rotatory  -inert -a  on  cannon 
vibrations,  three  uniform  cross  section  pipe  structures  resembling  contemporary 
tank  cannons  were  analyzed  using  this  technique.  The  resu^s  of  this  analysis 
are  contained  in  Figure  2  where  the  percent  error  (with  respect  to  the  full 
Timoshenko  equation)  is  plotted  against  mode  number  for  the  three  structures 
analyzed.  The  results  indicate  that  both  of  the  simpler  models  predict  fre¬ 
quency  values  which  are  higher  than  those  which  the  Timoshenko  equation  would 
yield.  The  'thicker'  structures  are  less  accurate  when  the  Euler-Bernoull i 
equation  is  used,  however,  the  inclusion  of  rotatory  inertia  in  the  Timoshenko 
equation  contributes  very  little  for  any  of  the  structures.  At  an  outer 
diameter  of  8  inches,  a  7.5  percent  error  in  frequency  would  occur  in  the  sixth 
mode  using  the  simplest  formulation,  while  only  1.5  percent  error  would  result 
using  the  complex  equation  without  rotatory  inertia. 

ERROR  ESTIMATES  IN  FREQUENCY  CALCULATIONS  FOR  UNIFORM  BEAMS 


PERCENT  ERROR  VS  MODE  NUMBER 


Figure  2.  Error  in  frequency  calculations  for  uniform  beams. 
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A  comparison  of  mode  shapes  for  the  sixth  mode  of  the  8-inch  beam  is 
plotted  in  Figure  3.  The  upper  curve  represents  the  shape  using  the  Timoshenko 
equation,  while  the  lower  is  the  same  for  the  Euler-Bernoul 1 i  equation.  As  can 
be  seen,  they  are  quite  similar  except  at  the  end  points  where  a  deflection 
error  of  -8.3  percent  and  a  slope  error  of  -0.75  percent  results  from  the  use  of 
the  simpler  equation.  Although  this  appears  to  be  quite  substantial,  these 
error  levels  become  less  for  the  lower  modes.  Since  mode  shape  contribution 
diminishes  with  increasing  mode  number,  these  error  estimates  become  relatively 
unimportant  to  the  overall  transient  response.  In  addition,  a  220-inch  long  by 
8-inch  diameter  105-mm  cannon  is  heftier  than  any  fielded  or  developmental  hard¬ 
ware,  therefore,  the  reported  errors  become  upper  bounds  when  actual  cannons  are 
modelled.  For  these  reasons,  the  Euler-Bernoull i  formulation  of  the  vibration 
equation  was  selected  for  the  remainder  of  this  analysis. 


MODE  SHAPE  COMPARISON  6TH  MODE  8"  0D  BEAM 
SHAPE  VS  LENGTH 


Figure  3.  Mode  shape  comparison  for  uniform  beam  models. 


Modelling  the  Non-Prismatic  Effect 


It  goes  without  saying  that  the  inclusion  of  non-prismat ic  geometries 
increases  the  complexity  of  this  problem.  The  normal  vibration  functions  for 
the  uniform  beam  were  shown  to  be  analytic  for  any  of  the  model  equations 
chosen.  This  is  a  great  advantage  in  solving  transient  problems  where  the  modal 
functions  need  to  be  differentiated  for  certain  types  of  driving  loads.  (This 
is  especially  true  for  gun  tube  loads  as  will  be  shown  later.)  If  the  modal 
function  is  approximated  by  a  set  of  tabular  values  (e.g.,  numeric  tables  from 
finite  element  analysis),  accurate  differentiation  is  difficult,  even  when  the 
points  are  connected  by  approximating  polynomials.  Since  the  non-uniform  nature 
of  a  gun  tube  is  consistent  in  that  it  is  composed  of  cylindrical  sections, 
tapered  sections,  and  step  changes  between  sections,  attempting  to  solve  the 
free  vibration  equation  with  a  single  analytical  function  may  be  impossible. 

A  number  of  approaches  were  found  during  the  research  phase  of  the  study 
(refs  9-16).  Each  of  these  authors  attempted  to  resolve  the  problem  of 
modelling  the  flexural  dynamics  of  a  non-prismatic  beam.  All,  however,  fell 
short  of  expectations  in  terms  of  fulfilling  the  need  for  analytically 
expressible  mode  shape  functions  differentiable  through  the  second  order. 

Canned  finite  element  routines  were  ruled  out  early  due  to  the  lack  of  ver¬ 
satility  in  their  ability  to  model  loads  which  are  functions  of  the  dependent 
variable.  In  regard  to  their  employment  for  determining  the  mode  shapes,  the 
only  way  a  continuous  analytical  function  could  be  achieved  was  through  a  least 
squares  or  spline  fit  using  the  normal  modes  displacement  information.  This  was 
unattractive  since  it  is  well  known  that  modal  displacements  approach  sinusoidal 
shapes  for  beams  with  mild  non-uniformities.  Any  method  developed  for  struc¬ 
tures  of  this  type  should  attempt  to  exploit  this  property. 


Upon  conversing  with  and  reviewing  the  work  of  Dr.  Sneck  (ref  17),  a  method 
of  dealing  with  the  problem's  non-prismatic  nature  was  conceived.  He  suggested 
that  the  tube  be  divided  into  two  uniform  segments  of  differing  cross  sections, 
each  of  which  is  analyzed  using  the  Euler-Bernoull i  equation.  By  the  separation 
of  variables  technique  (refs  2,3),  the  application  of  the  appropriate  boundary 
conditions,  and  the  enforcement  of  interfacial  continuity  at  segment  boundaries, 
expressions  involving  trigonometric  and  hyperbolic  terms  for  the  mode  shape 
functions  will  result.  Each  segment  will  have  its  own  set  of  mode  dependent 
coefficients.  Continuity  across  segment  boundaries  will  be  accomplished  by 
equating  boundary  values  of  displacement,  slope,  bending  moment,  and  shear  as 
calculated  by  each  adjacent  mode  function.  This  method  seemed  ideal  for  the 
problem  at  hand,  and  it  was  chosen  contingent  upon  its  adaptability  to  cannon 
models  where  more  than  two  segments  would  be  needed.  It  was  felt  that  any 
number  of  segments  could  be  handled  by  this  method,  hence,  the  term  "Uniform 
Segment  Method"  (USM)  was  coined. 

The  equation  defining  the  mode  shapes  for  a  general  case  employing  the  USM 
is 

j  (x )  =  A^j  cosh  a-jjx  +  B-jj  sinh  a-jjx  +  C-jj  cos  a-jjx  +  D-jj  sin  a-jjx  (10) 

where: 

i  =  segment  number 

j  =  mode  shape  number 

x  =  normalized  length  (x/L) 

A , B , C , D  =  mode  shape  coefficients 

r  1  w  . 

O;  ;  =  argument  coefficient  [Lru>;(-  (--)■;)*] 

J  J  g  El 

By  employing  this  segmented  mode  shape  function  and  the  well-known  separation  of 
variables  technique  on  the  pure  bending  equation,  the  unknown  coefficients  and 


frequency  for  each  mode  may  be  tracked.  For  a  model  utilizing  M  sections,  the 
imposition  of  the  free-free  boundary  conditions  will  yield  the  following 
equations  relating  bending  moment  and  shear  at  the  extremities  of  the  beam: 


I  I 


4>l,j(0) 

=  0 

(11a) 

*MJ(1> 

=  0 

(lid) 

♦ij(°> 

=  0 

(lie) 

=  0 

(lid) 

Continuity  at  the  interfacial  locations  is  preserved  by  equating  the  values  of 
displacement,  slope,  bending  moment,  and  shear  as  calculated  by  adjacent  mode 
functions  at  the  interfacial  boundaries  between  segments.  These  conditions  may 
be  written  as  follows: 


♦i-l,j<xi-l>  =  *i,j(*i-l>  (12a) 

♦i-l.jtXi-l)  *  *i,j(*i-i>  (12b) 

(El )  ( j  (*i-i )  =  (El),*- jfXi.j)  (12c) 

(EDt.jo’j  ,j(xi-l>  «  (El )  j(Xi_i )  ( 1 2d ) 


where: 

i  =  segment  number 

=  normalized  axial  coordinate  of  segment  k's  right  boundary 
For  a  model  containing  M  segments,  there  are  4M  coefficients  and  one  frequency 
to  be  evaluated  for  each  mode  shape.  By  setting  =  l  renders  the  system 

deterministic  for  which  unique  solutions  exist  at  each  natural  frequency.  A  set 
of  algebraic  equations  in  the  natural  frequency  (wj)  and  the  4M-1  unknown  mode 
shape  coefficients  is  presented  in  matrix  form  in  Figure  4.  The  system  matrix 
is  shown  to  contain  an  orderly  set  of  entries  representing  the  boundary  con¬ 
ditions  (B.C.)  and  matching  conditions  (M.C.).  A  2X4  subarray  beginning  at 


global  location  1,1  contains  the  terms  for  evaluating  the  B.C.  at  x  =  0,  whereas 
a  similar  array  ending  at  N,N  (N=4M)  represents  the  B.C.  at  x  =  1 .  'Walking1 
along  the  main  diagonal  are  M-l  subarrays  4X8  in  size  which  represent  the  pres¬ 
ervation  of  continuity  across  segment  boundaries.  The  remaining  terms  in  the 
matrix  are  zero.  The  non-zero  entries  are  functions  of  the  segment  properties, 
boundary  location,  and  the  unknown  natural  vibration  frequencies.  Setting  the 
determinant  equal  to  zero  and  solving  for  the  roots  of  the  resulting  charac¬ 
teristic  equation  produces  values  for  these  frequencies.  At  each  frequency,  a 
reduced  system  is  developed  by  setting  Ajj  =  l,  eliminating  the  first  row,  and 
shifting  the  first  column  to  the  right  side  of  the  equation.  The  solution  to 
this  linear  system  yields  the  remaining  normalized  mode  shape  coefficients. 


N-SEGMENT  UNIFORM  PROPERTIES  SYSTEM  OF  EQUATIONS 

Lil  (2  X  4) 

B.C.  «f  XrO 


M.C.  of  X  =  Xj 
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ZERO 


M.C.  of  X  =  X2 
(4  X  8) 


ZERO 


N  X  N  SYSTEM  MATRIX 


N  X  f 

COEFF  I  C  I  ENT  VECTOR 
M  =  M/4 

I  -  SEGMENT  NUMBER 
|  -  MODE  NUMBER 


Figure  4.  Uniform  segment  system  of  equations. 
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Unlike  finite  element  techniques,  this  method  aoes  rot  reo^-ire  the 
segmenting  of  uniform  sections.  Rather,  they  transform  cn  a  one-to-one  basis 
from  the  real  to  the  model  plane,  thus  only  a  few  segments  are  needed  to  accu¬ 
rately  represent  gently  varying  non-uniformities.  Additionally,  the  number  of 
modes  available  for  calculation  is  not  a  function  of  the  number  of  segments.  In 
finite  elements,  the  mode  shapes  are  defined  by  the  disolaced  locations  of  the 
model's  nodes.  There  are  not  enough  points  available  for  an  accurate  represen¬ 
tation  of  the  mode's  shapes  at  higher  natural  frequencies.  This  is  not  the  case 
for  the  USM,  since  nodeless  elements  using  trigonometric  and  hyperbolic  terms  in 
the  mode  shape  functions  are  employed. 

To  determine  the  accuracy  of  this  method  for  free  vibration  modelling  of 
gun  structures,  the  predictions  from  an  established  finite  element  code  (ABAQUS 
(ref  18))  for  a  typical  gun  tube  (105-mm  M68)  were  compared  to  the  results 
generated  by  the  dedicated  routines  (MODE:)  written  by  the  atj'hor  in  support  of 
this  analysis.  The  gun  and  its  modelling  representations  are  shown  in  Figure  5. 

The  M68  gun  tube  is  210  inches  long  containing  two  tapered  and  three 
cylindrical  sections.  The  physical  schematic  is  shown  in  the  top  sketch  in 
Figure  5.  This  finite  element  model  (FEM)  was  comprised  of  54  equally-spaced 
nodes  with  two  degrees  of  freedom  (dof)  per  node  and  53  prismatic  beam  elements. 
The  sketch  in  the  middle  of  the  figure  indicates  the  element  density  in  each 
physical  section  of  the  structure.  On  the  lower  sketch  of  the  same  figure, 
segmentation  of  the  USM  model  is  shown.  Only  six  segments  of  four  dof  each  were 
used  for  a  total  of  24  dof.  The  natural  frequencies  and  mode  shapes  for  the 
first  six  vibration  modes  were  calculated  using  both  the  Timoshenko  and  Euler 
equations  for  the  FEM.  In  the  USM  analysis,  only  the  Euler  equation  is 


employed. 


1 05mm  M 6 8  GUN:  COMPARISON  OF  ANALYTICAL  MODELS 
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CUN  SCHEMATIC:  210  INCHES  LONG  (NUMBER  OF  INCHES) 


FINITE  ELEMENT  MODEL:  108  D  0  F  (NUMBER  OF  ELEMENTS) 


UNIFORM  SEGMENT  MODEL:  24  D  0  F  (NUMBER  OF  SEGMENTS) 


Figure  5.  105-mm  M68  Gun.-  Comparison  of  analytical  models. 

The  results  of  the  analysis  are  shown  in  Figures  6  and  7.  In  Figure  6  a 


comparison  among  the  frequency  calculations  is  presented  in  the  form  of  a  bar 


graph.  The  abscissa  contains  the  first  six  mode  numbers,  while  the  ordinate  is 


labelled  with  the  Log^o  value  of  the  frequency.  As  is  indicated,  the  FEM  using 
the  Timoshenko  equation  produces  the  lowest  values  for  all  six  frequencies, 
while  the  USM  produces  the  highest.  The  differences  in  the  extremes  are  between 
six  and  ten  percent  with  the  greatest  discrepancy  occurring  at  the  highest  mode 
number . 
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105mm  M68  GUN:  NATURAL  FREQUENCY  ESTIMATES 


Figure  6.  105-mm  M68  Gun:  Natural  frequency  estimates. 


105mm  M68  GUN:  MODE  SHAPE  COMPARISONS 


NORMALIZED  SHAPE  OF  6TH  MODE 


Figure  7.  105-mm  M68  Gun:  Mode  shape  comparison. 
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In  Figure  7  the  normalized  mode  shape  plots  for  the  sixth  mode,  as  calcu¬ 
lated  by  each  of  the  models,  is  shown.  The  sinusoidal  nature  of  the  mode  shape 


To  implement  these  effects,  an  independent  set  of  routines  ( I N I T : )  was 
developed  to  generate  the  initial  conditions  for  the  gun/beam.  These  routines 
accept  the  geometry  from  the  normal  modes  analysis,  adding  two  resilient  sup¬ 
ports,  point  loads  simulating  non-structural  mass,  and  distributed  loads  repre¬ 
senting  each  segment's  weight.  The  deflection  response  of  the  supports  may 
include  clearance  and  non-linear  elasticity.  Initially,  the  rigid  body  solution 
for  the  gun/beam  is  solved.  This  yields  support  reactions,  deflections,  and  the 


if 


where 


=  i-th  mode  shape 


=  total  number  of  modes 


=  k-th  uniform  segment 
=  total  number  of  segments 
=  m-th  applied  load 
s  total  applied  loads 
=  bending  resistance  of  k-th  segment 
=  weight/length  of  k-th  segment 
=  i-th  mode  shape  of  k-th  segment 


q-j  =  i-th  modal  amplitude 

pm(x,t,...)  =  m-th  applied  load  function 


(EDk 


wk(y) 


=  static  mass  load  applied  to  k-th  segment 


=  static  deflection  curve 


Since  the  modal  amplitudes  are  time-dependent  only  and  the  static  deflec¬ 
tion  is  space-dependent  only,  the  above  equation  may  be  separated  as  follows: 


\  wi(  r  r 

<EI)k  (  ^  ♦kifli)  +  g-  (  }  4>ki<*i)  =  /  Pm(x. 


(EDk  ys  «  Wk(y)  (H 

The  solution  to  Eq.  (14a)  yields  the  dynamic  response  with  respect  to  the 
initial  shape  of  the  beam.  The  solution  to  the  second  equation  is  the  static 
deflection.  Superposition  of  these  leads  to  the  total  dynamic  response  with 
respect  to  a  global  inertial  coordinate  system.  Since  Eq.  (14b)  has  already 
been  solved  (last  section),  the  solution  to  Eq.  (14a)  will  close  the  problem. 


By  the  nature  of  the  segment  mode  shape  functions 
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4*ki  =  (-[')  4>ki 


which  leads  to 


5  ♦ki 


i  =  l 


.El.  .aki,4 

Pi  +  9(--)k(-[-)  Pi 


g_ 

*k 


5 

m=l 


Pm(*»t, • ■ • > 


Recalling  a  parameter  from  Eq.  (10)  and  employing  current  notation 


(15) 


(16) 


«ki  -  u-«,  (|lgI3k)*‘ 


which  upon  substitution  yields 

I 


/  ♦ki [Pi  +  Wi2Pi3  =  J- 
i  =  l 


^  pm ( x , t , . . . ) 

m=l 


(17) 


For  non-dissipative  boundary  conditions,  the  orthogonality  of  the  mode 
shapes  with  respect  to  the  weight  function  w(x)/g  leads  to  the  following: 


L  ~g~^  «Mj  =  }  /  g"  ♦ki^kj  dX  = 

0  9  J  *-  x^o  9 

k=l  KU 

0  for  i  *  j 
Gj  for  i  =  j 

where 

xk0  =  axial  location  of  lower  boundary  of  segment  k 
X)<2  =  axial  location  of  upper  boundary  of  segment  k 
Note:  xkl  =  x(k+1)0 
w(x)  =  structural  weight  distribution 


(18a) 

(18b) 
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By  multiplying  Eq.  (17)  by  ----  integrating  over  the  length  and 
applying  Eq.  (18),  the  following  ordinary  differential  for  modal  amplitudes 
results : 

M  K 

1  v  v  .xkl 

9i  +  Ui2Pi  =  ?-  )  )  /  pn,(x,t, . . .  )<J)ki  dx  (19) 

1  1  L  I  Xk0 

m=l  k=l 

The  solution  to  this  equation  yields  a  vector  of  modal  amplitudes  used  to 
calculate  the  dynamic  displacement,  velocities,  and  slopes.  These  equations  are 
as  follows: 


y<j(x)  =  }  q i  ( t  )4>i  ( X ) 


(20a) 


yd(x)  =  ^  q-j(t)<Mx) 


(20b) 


y'  (x)  =  ])  q-j  ( t )<t>i  ’  (x) 


(20c) 


where 


-i  ( x )  =  4>ki(x^)  segment  mode  shape 


Representation  of  the  Forcing  Functions 

Previous  analyses  and  test  results  indicate  that  static  droop  dominates  the 
overall  curvature.  Dynamic  response  is  shown  to  be  a  displacement  perturbation 
about  the  static  shape.  The  maximum  displacement  levels  are  an  order  of  magni¬ 
tude  less  than  the  static  muzzle  deflection.  Since  the  loading  functions  unique 
to  gun  dynamics  are  dependent  upon  a  gun’s  overall  shape,  so  is  the  dynamic 


response.  Recalling  the  functional  form  for  these  load  functions 


PmfX't-y-y’ -y-y-y".y' 1 


from  Eq.  (14a),  the  following  expression  from  Ec.  (19) 
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•  •  )4>ki  dx 


must  be  evaluated  for  each  load  function  considered.  These  load  functions  and 
evaluations  will  now  be  presented  and  discussed. 

1.  Recoil  inertial  load.  An  inertia  couple  develops  within  any  segment  k 
during  recoil  which  is  expressible  as 

Wl- 

Pl ( t . x , y ’ )  =  --  ar(t) ( (l-x)y’ ] ’  (21 

where 

ar(t)  =  recoil  acceleration 
i  =  total  length  of  gun 

y  =  total  transverse  displacement  (static  and  dynamic) 

Details  of  the  derivation  may  be  found  in  Reference  20.  Evaluation  of  the 
expression 

1  A  wk 

t  /  --  ar  ( t )  t  (l-x)y’  3  '<J>ki  dx 

o  9 

upon  the  substitution  of 

I 

y  =  ^  4>kjdj  +  ys  (22 

j=l 

and  normalizing  with  respect  to  tube  length  yields  the  following: 


where 


I 

1,1  ,  r  1  '  xkl 

q~  /  Pl(*.t,y'  )4>ki  dx  =  Ki(t){wk^kj[(l-x)(  2  <J>kjdj+ys)]- 

j=l  *kj 
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(23) 


x  =  x/J 

2.  Pressure  curvature  load.  Due  to  the  nature  of  curvature  within  real 
bodies,  diametrically  opposite  bore  'surfaces'  possess  differing  'areas'.  A 
pressure  load  acting  within  this  confined  chamber  will  tend  to  straighten  the 
chamber  because  the  area  of  the  concave  'surface'  is  greater  than  its  convex 
counterpart.  The  expression  for  this  load  function  (ref  20)  is 

P2(x,t,y")  =  -ABPB(t)y"{H(xp(t)  -  x) )  (24) 

where 

Ab  =  bore  area  of  tube 
PB(t)  =  propellant  gas  pressure 
Xp(t)  =  projectile  location 
H(Xp-x)  =  step  function 

Upon  substitution  from  Eq.  (22),  the  evaluation  of  the  normalized  integral 
yields 


P2(x.t,y")<J)ki 


dx 


>1 


S 


(25) 


3.  Projectile  trajectory  load.  The  accelerating  projectile,  although  of 
considerably  less  mass  than  the  tube,  can  exert  a  significant  transverse  force 
when  it  is  confined  to  travel  along  a  curved  path.  Simkins  (ref  20)  identified 
this  loading  and  derived  the  following  expression.- 

w 

P3(x,t,y' ,y",y)  =  -  [y  +  2xp(t)y'  ♦  (xp(t))*y"  +  g] 

6(xp(t)-x)  (26 

where 

wp  =  projectile  weight/unit  length 
6(xp(t)-x)  =  Dirac  delta  function 

Upon  substitution  of  Eq.  (22)  and  normalizing  with  respect:  to  tube  length, 
the  integral  expression  becomes 

1. 1  •  -  - 

p-  /  P3(*.t,y'  ,y",y)<J>ki  6(xp(t)-x)dx  = 

ui  o  y 
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(*P)M  2  ^j(xp^j  +  ys(xp)l  +  g)  (27 

j=i 


where 


Wp  =  projectile  weight 


*J 

l] 
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4.  Projectile  eccentricity.  Wu  (ref  23)  has  postulated  the  existence  of  a 
couple  exerted  on  the  tube  when  the  projectile  propelling  force  does  not  pass 
through  the  projectile's  mass  center  which  is  expressible  as 


P4<x,t)  =  ABFB<t)(  (Xp-x)sin(e0+2rrrxp) 


where 


ep  =  radial  eccentricity  of  projectile 
Jp  =  wheelbase  of  projectile 
r  *  rifling  twist  (=  0  for  smooth  bore) 


/  6' (xD-x)dx  =  1 


Upon  substitution  of  Eq.  (22)  and  normalizing,  the  integral  expression 


becomes 
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r  J  P4  (x ,  t  )4>ki  dx  =  -K4  ( t )  4>ki(xp)  (29] 

ui  o  1 

where 

ep 

K4(t)  =  (^-)sin(e0+2rtTxp)i 

°i  *p  K 

5.  Projectile  rotational  load.  As  the  projectile  travels  along  a  rifled 
tube,  its  mass  center  rotates  causing  a  centrifugal  load  on  the  bore.  The  ver¬ 
tical  component  of  this  reaction  is  expressible  as 


P5 (x , t ) 


[(27TTXp)z6(Xp-X)  sin(0o  +  27TTXp) 


(2nTxp)cos(4>0+2nTXp)  ]6(xp-x) 


where 


=  projectile  axial  velocity 
*  projectile  axial  acceleration 
=  rifling  twist  in  revolutions  per  inch 


Upon  substitution  and  normalization,  as  above,  ;ne  Integra'*  evaluation 
becomes 

1,1 

z-  I  p5(x,t)0ki  dx  =  K5  ( t  )<t>j,  .j  (xD) 

“to  h 

where 

K5  ( t  >  =  [  ( 27TTXp )  2si n ( 0o+2rtTXp)  -  ( 2rxp  Icos  {4>0+2my.-, )  ] 

6.  Stationary  mass  reaction.  The  stationary  non-structural  masses 
(breech,  bore  evacuator,  muzzle  brake)  cause  transverse  inertial  loadings 
expressible  as 


where 


P6(*'t<y)  =  -  --  y6fxs-x) 


ws  =  weight  of  stationary  mass/length  of  mass 
xs  *  axial  location  of  mass 

Xs  =  xs/ 

The  integral  evaluation  is 


where 
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Ws  =  weight  of  stationary  mass 

7.  Stationary  mass  eccentricity.  Should  these  masses  possess  vertical 
eccentricity  with  respect  to  the  bore  axis,  a  couple  will  develop  causing  a 
expressible  as 


(31) 


(32) 


(33) 


load 


P7<x-t)  *  -  ABPB ( t ) ( "  ) ^ ' (xs~x) 

"R  *s 


(34) 


where 


Wr  =  total  recoiling  weight 

es  =  vertical  eccentricity  from  bore  axis  (  +  =>  above  axis) 


ls  =  wheelbase  of  eccentric  mass 

The  integral  evaluation  of  this  load  yields 


1  J  • 

5-  P7(x.*)4>ki  dx  =  K74>k-i  < xs ) 


o 


(35) 


where 


K7  =  -  g'Wr  Wt)(es> 


8.  Support  reactions.  The  mounting  supports  may  be  characterized  as  non¬ 
linear  spring  elements  which  react  actively  with  the  total  displacement  of  the 
tube.  Initially,  the  support  deflection  balances  the  static  loads.  As  the  tube 
vibrates,  the  displacements  of  the  tube's  support  locations  change.  For  a 
general  non-linear  support  spring,  the  reactive  load  may  be  expressed  as 

p8(x,t,y)  =  f>n(y>6(Xrn-x)  (36) 

where 

Frn(y)  =  dynamic  reaction  force  of  n-th  support  (total  reaction  -  static 
reaction ) 

xrn  =  axial  location  of  n-th  support 
The  integral  equation  is  as  follows: 


1  J 

5:  /_  Pe(x,t,y)<t>ki  dx  =  Kg(y)$ki  (xrn) 


0 


(37) 


where 


k  ,  ,  Fm(V) 
k8  ( y )  =  ~g;~ 
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9.  Bore  eccentricity  load.  Should  the  gun's  bc'c  ar c  oute-  c  aneter  be 
non-concentric ,  each  differential  segment  of  the  bean  contributes  an  inert'* 
couple  about  the  axis  which  may  be  distributed  in  va-ying  Degrees  over  the  full 
length  of  the  gun.  The  differential  load  per  unit  length  may  be  expressed  as 


fol lows : 


where 


Pg(x,t)  =  -wkar(t)e(xe)6' (xe-x ) 


ar(t)  =  o~  :  normalized  acceleration 

r  WR 


e(xe)  =  transverse  eccentricity  at  xe 

Since  this  is  a  differential  load,  its  total  contribution  becomes 

1  ,1  .... 
r~  J  Pg(x,t)$k.jdx  =  -Kg(t)e(xe)<t>.,  (xelAx 

ui  0 


where 


wkar ( t ) 1 

........ 


Ax  =  j-  =  differential  length  over  which  the  eccentricity  is  distributed 
The  total  load  is  the  summation  of  the  differential  loads.  When  the  formulation 
is  cast  in  a  continuous  form,  the  final  value  for  the  load  becomes 

1  *  .e(x)$i(x).  X1  *1  -  - 

/  Pg  ( x ,  t  )<t>k-jdx  =  -Kg  ( t )  { - | - 1-  -  /-  e'  (x)4>,  (x)dxj  (40) 

i  0  1  *0  x0 

Hodal  Amplitude  O.O.E.'s  and  Numerical  Solution  Process 

The  algebraic  rearrangement  of  the  loading  functions  developed  in  the  pre¬ 
vious  section  results  in  a  system  of  O.O.E.'s  in  the  amplitude  vector  q ^ ( t ) . 

The  general  equation  is 


(t )  ♦  [C]q-j  ( t )  ♦  [K]qi(t)  = 


/  f mi ( * ) 


where 


[M]  =  inertia  matrix  of  order  I 
[C]  =  Coriolis  matrix  of  order  I 
[K]  =  stiffness  matrix  of  order  I 
fmi(t)  =  i-th  mode;  m-th  load  driving  force 

The  three  matrices  are  fully  populated  unlike  the  case  of  a  discrete 
spring,  mass,  and  damper  system  in  which  the  inertia  matrix  is  diagonal.  For 
the  case  concerning  gun/beam  dynamics,  the  projectile  and  stationary  non- 
structural  masses  cause  coupling  between  mode  shapes  showing  up  as  off-diagonal 
terms  in  the  inertia  matrix.  The  Coriolis  matrix  is  appropriately  named  because 
its  only  contributing  load  is  due  to  the  moving  projectile  travelling  along  a 
moving  path  created  by  the  vibrating  tube.  A  Coriolis  force  component  results 
from  this  interaction.  Modal  coupling  is  a  characteristic  of  this  load,  there¬ 
fore,  matrix  [C]  is  fully  populated.  The  stiffness  matrix  [K]  contains  the 
natural  vibration  frequencies  (w^)  along  the  main  diagonal  as  well  as  other 
contributions  from  recoil  inertia,  pressure  curvature,  and  projectile  trajectory 
forces.  Inclusion  of  these  terms  causes  the  matrix  to  be  fully  populated. 

The  forcing  functions  on  the  right  side  of  the  equation  are  all  time- 

dependent  in  that  the  ballistic  force,  recoil  inertia,  or  projectile  location 
and  kinematic  state  are  needed  for  their  evaluation.  The  static  slope  is 
required  for  recoil  inertia  load  evaluation,  whereas  curvature  is  needed  in  the 

pressure  curvature  and  projectile  trajectory  loads.  Point  loads  such  as  those 

due  to  projectile  location  travel  along  the  structure,  while  those  due  to  the 
stationary  masses  are  fixed  in  the  spatial  coordinate.  The  support  reactions 
require  an  evaluation  of  the  total  deflection  of  the  tube  at  their  fixed  loca¬ 
tions  with  the  reaction  force  being  a  function  of  the  total  displacement. 


From  a  solution  standpoint,  the  problem  is  qu'te  cox^pv  due  to  the 

inertia  coupling.  In  order  to  solve  the  system  of  eauaticns  using  numerical 
procedures,  the  matrix  [M]  must  be  triangularizec  r y  an  elimination  process 
(Gauss)  with  back  substitution  performed  to  arrive  at  a  solution  to  the  modal 
acceleration  vector  q^(t).  Rearrangement  of  Eq.  (41)  and  introduction  of  a 
discrete  time  step  yields 

M 

[M]di  (tn)  =  J  fmi(tn)  -  [C]q,(tn)  -  [k]q -j  ( tn )  (42) 

m=l 

where  tn  is  any  integration  time.  Initially, 

t0  =  0  and  q-j  (t0)  =  q ( t0 )  =  0  (43) 

These  conditions  allow  for  the  startup  of  the  solution  process.  Equation 
(42)  is  solved  for  q^(tn)  by  back  substitution  into  the  tr iangularized  inertia 
matrix.  The  resultant  acceleration  vector  is  integrated  to  yield  the  velocity, 
and  finally  the  modal  displacement  vectors. 

A  predictor-corrector  technique,  based  upon  the  Adams-Bashf orth-Moul ton 
multi-step  formulation  (ref  24)  adapted  for  systems  of  equations,  was  the 
algorithm  chosen  for  the  integration  process.  A  fixed  time  step  with  con¬ 
vergence  control  and  limited  iteration  steps  is  provided  by  the  user  through 
computer  input  files.  This  multi-step  method  needs  four  starting  values  of  the 
function  being  integrated.  These  values  are  generated  by  using  a  Taylor  series 
approximation  to  the  solution  equation  through  the  four  initial  time  steps.  The 
predictor  portion  provides  an  initial  solution  for  the  modal  velocity  and 
displacement  amplitudes  by  using  an  integrating  algorithm  based  upon  the 
Adams-Bashforth  Four-Step  Method.  This  is  an  explicit  method  requiring  function 
evaluations  from  four  preceding  time  steps.  The  solution  predicted  is  used  in 
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the  corrector  portion  of  the  algorithm  which  is  an  implicit  technique  known  as 
the  Adams-Moul ton  Three-Step  Method.  Iterations  on  the  approximate  solution 
vector  continue  until  convergence  is  assured  based  upon  a  criterion  supplied  by 
the  user.  If  the  criterion  cannot  be  met  in  the  maximum  number  of  iteration 
steps,  the  computer  routines  report  this  occurrence  allowing  for  user  interven¬ 
tion.  Additionally,  an  error  estimate  including  the  number  of  iterations 
attempted  is  available  as  output  for  each  integration  step.  This  is  useful  for 
determining  a  suitable  time  step  for  a  given  analysis. 

RESULTS  AND  DISCUSSION 

At  this  time,  these  modelling  routines  are  being  debugged  and  tested 
against  the  predictions  from  standard  solutions  and  other  independent  analyses. 
Subsequently,  this  modelling  will  be  compared  with  experimental  data  (both  field 
and  laboratory  generated)  to  identify  any  inherent  shortcomings.  Parametric 
studies  addressing  the  projectile's  exit  vector  (see  Figure  7,  Reference  1)  and 
its  sensitivity  to  perturbations  in  the  design  and  operational  parameters  of 
fielded  weapons  will  follow. 
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